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Abstract 

Eg x E$ heterotic string and M-theory, when compactified on smooth 
Calabi-Yau manifolds with 577(4) vector bundles, can give rise to 
softly broken N = 1 supersymmetric theories with the exact matter 
spectrum of the MSSM, including three right-handed neutrinos and 
^ one Higgs-Higgs conjugate pair of supermultiplets. These vacua have 

O the SU(3)c x SU(2) L x U(l)y gauge group of the standard model 

augmented by an additional gauged U(1)b-l- Their minimal content 
• • requires that the B-L symmetry be spontaneously broken by a vacuum 

. ^ expectation value of at least one right-handed sneutrino. The soft su- 

^ persymmetry breaking operators can induce radiative breaking of the 

B-L gauge symmetry with an acceptable i?-L/electroweak hierarchy. 
In this paper, it is shown that U(1)b-l cosmic strings occur in this 
context, potentially with both bosonic and fermionic superconductiv- 
ity. We present a numerical analysis that demonstrates that boson 
condensates can, in principle, form for theories of this type. However, 
the weak Yukawa and gauge couplings of the right-handed sneutrino 
suggests that bosonic superconductivity will not occur in the simplest 
vacua in this context. The electroweak phase transition also disallows 
fermion superconductivity, although substantial bound state fermion 
currents can exist. 
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1 Introduction 



Heterotic string and M-theory when compactified on smooth geometric and 
vector bundle backgrounds [I]-[H] can give rise to "heterotic standard mod- 
els" PHI dH H2]; that is, four-dimensional N = 1 supersymmetric theories 
with exactly the matter and Higgs spectrum of the MSSM. Supersymmetry 
can be spontaneously broken by non-perturbative effects in the hidden sector. 
Integrating out this sector, the low-energy theory contains "soft" supersym- 
metry breaking operators whose generic form is well-known p2]-[T8]. To be 
phenomenologically viable, any such theory must have two properties: 1) 
three right-handed neutrino chiral multiplets, one per family, and 2) "matter 
parity", a discrete Z2 symmetry which prohibits too rapid baryon and lepton 
number violating processes [T9]-[22]. 

These two properties are most easily satified in heterotic standard models 
constructed using vector bundles with SU (4) structure group [23] -[2E]- In ad- 
dition to the MSSM spectrum, such vacua have three right-handed neutrino 
chiral multiplets, thus satisfying the first property. The low-energy gauge 
group also contains the SU(3)c x SU(2)l x Z7(l)y gauge group of the stan- 
dard model augmented, however, by a gauged U(1)b-l factor. This contains 
matter parity as a discrete subgroup. If the B-L symmetry could be spon- 
taneously broken to its matter parity subgroup, then the second property 
would be satisfied as well. However, this is not possible in smooth heterotic 
compactifications since the necessary 3(B — L)-even multiplets are disallowed 
as zero-modes. A second solution is to have U(1)b-l radiatively broken at 
low energy, not too far above the electroweak scale. It would then act as a 
custodial symmetry for matter parity, suppressing baryon and lepton violat- 
ing decays yet not unduly affecting electroweak physics. 

In several recent papers [301 EI], it was shown using a quasi-analytic 
solution to the renormalization group equations (RGEs) that this can indeed 
occur for a range of initial soft breaking parameters. Scaling down from 
the compactifcation mass, the gauged U{1)b-l is first spontaneously broken 
by a non-zero vacuum expectation value (VEV) of the third family right- 
handed sneutrino. This is followed by radiative VEVs developing in the Higgs 
fields which induce an electroweak phase transition. The 5-L/electroweak 
hierarchy was found to be of C(10)-C(100). Recently, these results have been 
expanded to a much wider range of initial soft parameters using a completely 
numerical solution of the RGEs. This work will appear elsewhere [32] • Here, 
we simply note that this expanded range of parameters leads to three distinct 



possibilities for the soft mass squared parameters of squarks/sleptons at the 
electroweak scale. In addition to the negative third family right-handed 
sneutrino mass, 1) all such parameters are positive, 2) all are positive with 
the exception of a right-handed charged slepton and 3) all are positive with 
the exception of a left-handed squark. Each possibility can play an interesting 
role in cosmology 

The starting point of this paper is the assumption that smooth heterotic 
compactifications with SU(A) structure group are potentially phenomenolog- 
ically viable theories for low-energy particle physics. The distinct signature 
of this type of vacuum is that a gauged U{1)b-l symmetry is spontaneously 
broken at a low scale, not too far above the electroweak phase transition. As 
is well-known (33j El] , the breaking of a gauged Abelian symetry can lead to 
topologically stable cosmic strings. In principle, these strings can exhibit a 
wide variety of observable cosmological phenomena [55]- [60]. However, much 
of the analysis of cosmic strings has been carried out within the context of 
grand unified theories or specially constructed supersymmetric models whose 
spectra contain fields in addition to those of the MSSM with right-handed 
neutrinos [35]-[l3]. Furthermore, the coupling parameters associated with 
these fields are not constrained and can be assumed to be sufficiently large. 
As a rule, it is these extra fields that induce the potentially observable phe- 
nomena, such as bosonic or fermionic superconductivity |H] . In smooth 
B-L MSSM heterotic compactifications, there are no additional fields. The 
breaking of U{1)b-l and electroweak symmetry is accomplished via radiative 
expectation values for a right-handed sneutrino and Higgs fields respectively. 
As a consequence, the relevant parameters in this theory are those of the 
MSSM and, hence, tightly constrained by phenomenology. It follows that 
the existence and properties of cosmic strings in this context are severely 
restricted. 

In this paper, we analyze cosmic strings in the B-L MSSM theory. We 
show that such strings can indeed exist but are restricted to be BPS solu- 
tions at the critical boundary between Type I and Type II superconductors. 
There is a stable minimum of the scalar potential in which 1) B-L is broken 
by a VEV (^3) of the third family right-handed sneutrino and 2) electroweak 
symmetry is broken by Higgs expectation values (H°), (H°). At this min- 
imum, all other scalar fields have positive squared masses. Some, however, 
specifically right-handed charged sleptons and left-chiral squarks, have effec- 
tive masses (m 2 ) = m 2 + c(z/ 3 ) 2 with positive coefficient c, where m 2 is the 
associated soft breaking mass parameter. Three possibilites then arise. First, 
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all m? parameters might be positive, making the potential energy at the ori- 
gin of field space a minimum in all but the third family sneutrino direction. 
Second, a charged right-handed slepton could have a negative m 2 parameter, 
in addition to the right-handed sneutrino. Although the 5-L/electroweak 
vacuum remains the minimum, this destabilizes the slepton direction in the 
core of the cosmic string, potentially leading to a charge breaking conden- 
sate and bosonic superconductivity [UJ . The third possibility is that the soft 
mass squared parameter for a squark becomes negative, potentially leading 
to a charge and color breaking superconducting condensate. We have shown 
in [32] that each of these types of vacua is possible for a given range of initial 
soft parameters. From both the cosmological and phenomenological point of 
view, it is of interest to see if bosonic supercondictivity can occur in B-L 
MSSM cosmic strings. To explore this, we study a generic class of theories 
that arise in this context. Using a numerical analysis, a bound is derived 
that must be satisfied to allow the formation of a non-zero condensate and, 
hence, bosonic superconductivity. This analysis is then applied to the most 
straightforward B-L MSSM vacua and cosmic strings using simplifying as- 
sumptions. We find that the right-handed sneutrino Yukawa parameter and 
the Qb-l gauge coupling are too small to permit this essential constraint 
to be satisfied. We conclude that at least the simplest B-L MSSM cosmic 
strings do not exhibit bosonic superconductivity. 

B-L MSSM cosmic strings may also exhibit superconductivity induced by 
fermionic zero-modes in the string core [I5j SHI H7] . The fact that the gauged 
U {I) b-l extension of the MSSM is rendered anomaly free by the inclusion of 
three families of right-handed neutrino chiral multiplets plays an important 
role here. The cosmic string initially develops as a non-zero ra-fold winding of 
(b> 3 ) around some line in space. This couples directly to the left-chiral tauon 
ipE- and the chargino ipn+, forming left-moving fermion zero-modes on the 
string worldsheet. However, anomaly cancellation requires the appearance of 
right-moving fermionic modes, whose identity in the B-L MSSM context is 
not self-evident. We show that the third family left-handed sneutrino devel- 
ops a small VEV (JV3) following the electroweak phase transition. This wraps 
the core of the cosmic string with winding — n, opposite that that of (1/3). 
The right-chiral tauon ip e + and chargino ipH- couple to this field, inducing 
right-handed zero-modes which cancel all worldsheet anomalies. Thus, there 
is potential fermionic superconductivity in the cosmic string. We conclude, 
however, that the electroweak phase transition will, in general, render these 
fermionic currents unobservable |48j-|50j. 
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The paper is structured as follows. In Section 2, we briefly review the 
spectrum, superpotential and potential energy of the softly broken B-L 
MSSM theory. The structure of the B-L and electroweak breaking vac- 
uum is then presented in Section 3, including the effective scalar masses at 
this minimum and the 5-L/electroweak hierarchy. Section 4 is devoted to 
showing that the winding of the B-L charged right-handed sneutrino VEV 
around the origin leads to a cosmic string with critical coupling. The allowed 
patterns of soft scalar masses at the core of the cosmic string are discussed 
in Section 5. For each case, the stability criterion for a scalar condensate 
to develop in the string core, and, hence, for the string to be potentially 
superconducting, is derived. These criteria are then analyzed using a nu- 
merical analysis presented in the Appendix. Finally, in Section 6 we discuss 
potential fermionic zero-modes, show how the anomaly freedom of the B-L 
MSSM theory leads to appropriate left- and right-moving modes and present 
the constraints imposed on these currents by the breaking of B-L via the 
right-handed sneutrino. 

2 The N = 1 Supersymmetric Theory 

We consider an N = 1 supersymmetric theory with gauge group 

G = SU(3) C x SU(2) L x U(1) Y x U{1)b-l (1) 

and the associated vector superfields. The gauge parameters are denoted g%, 
g 2 , gy and qb-l respectively. The matter spectrum consists of three families 
of quark/lepton chiral superfields, each family with a right-handed neutrino. 
They transform under the gauge group in the standard manner as 

Qi = (3, 2, 1/3, 1/3), Ul = (3, 1, -4/3, -1/3), d % = (3, 1, 2/3, -1/3) (2) 

for the left and right-handed quarks and 

Li = (1,2, -1,-1), ^ = (1,1,0,1), ^ = (1,1,2,1) (3) 

for the left and right-handed leptons, where i — 1,2, 3. In addition, the spec- 
trum has one pair of Higgs-Higgs conjugate chiral superfields transforming 

as 

H= (1,2,1,0), # = (1,2,-1,0). (4) 
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When necessary, the left-handed SU{2) L doublets will be written as 

Qi = (Ui,Di), Li = (Ni, Ei), H=(H + ,H°), H = {H ,H~). (5) 

There are no other fields in the spectrum. The three right-handed neutrino 
chiral multiplets render this U(1)b-l extension of the MSSM anomaly free. 

The supersymmetric potential energy is given by the sum over the mod- 
ulus squared of the F and .D-terms. The F-terms are determined from the 
superpotential 

3 

W = fiHH + ^2 (K,iQiHui + Xd,iQiHdi + \ V)i LiEvi + X^UHei) , (6) 
t=i 

where we assume a mass-diagonal basis for simplicity. An innocuous mixing 
term of the form LiH as well as the dangerous lepton and baryon number 
violating interactions 

LiLjek, LiQjdk, Uidjdk (7) 

are disallowed by the U(1)b-l gauge symmetry. The SU(3)c and SU(2) L 
D-terms are of standard form. The U(l)y and U(1)b-l -D-terms are 

D Y = g Y <f ) \(Y/2) AB <PB (8) 

and 

D B -l = 9b-l4> ] a (Yb-l)^b (9) 

respectively, where index A runs over all scalar fields <pA- In the D-eliminated 
formalism, any Fayet-Iliopoulos parameters can be consistently absorbed into 
the definition of the soft supersymmetry breaking scalar masses. Hence, they 
do not appear in ^ and 

In addition to supersymmetric interactions, the potential energy also con- 
tains explicit soft supersymmetry violating terms. This breaking can arise in 
either F-terms, D-terms or both in the hidden sector. We will restrict our 
discussion to soft supersymmetry breaking scalar interactions arising exclu- 
sively from F-terms. Their form is well-known and, in the present context, 
given by [T5]-[T8"] 

^soft = ^ + ^3 S , (10) 
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where Vo, are scalar mass terms 



' 2* 
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V 2s = ^2(m 2 Q .\Qi\ 2 + m 2 Ui \ui\ 2 + m\\dif + m 2 Li \L t \ 2 + m 2 ^ 2 
i=i 

+m\\ei\ 2 + m^|#| 2 + m\\H\ 2 ) - (BHH + /ic) (11) 



and Vzs are the scalar cubic couplings 

3 

Vs s = ^(4*Qi#u< + A di QiHbi + A Vi LiHui + A £i L l He l + he). (12) 



i=l 



We choose the parameters in (11) and (12) to be flavor-diagonal. 



3 The B-L/ Elect roweak Hierarchy 

In [301 EI] a detailed one-loop renormalization group analysis of this theory- 
was carried out. In that analysis, tanfi was limited to 6.32 < tan/3 < 40 
and a specific range of initial parameters near the gauge unification scale 
M u ~ 3 x lO^GeK was chosen so as to allow a quasi-analytic solution of 
the RGEs. Here, we simply present the results. Subject to realistic, but 
constrainted, assumptions about the soft breaking parameters, it was shown 
that a hierarchy of radiative symmetry breaking takes place. 

First, at an energy scale of ~ TeV the third family right-handed sneu- 
trino soft mass parameter is negative; that is, m 2 3 < 0. It follows that this 
sneutrino acquires a vacuum expectation value (VEV) 

<-3> 2 = • (13) 

ilB-L 

Furthermore, evaluated at (z/ 3 ) all other scalars, including the Higgs fields, 
have vanishing VEVs. Therefore, this is a stable vacuum which sponta- 
neously breaks U{1)b-l while preserving the remaining SU(3)c x SU(2)l x 
U{1)y gauge symmetry. The Higgs effect then leads to identical masses for 
the B-L vector boson and the radial real scalar <5z/ 3 given by 

M Ab _ l = m Sua = V2g B _ L (u 3 ) . (14) 
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It is of interest to recall the expressions for the slepton/squark masses at this 
minimum of the potential. They were found to be 

I 2 \ 2 2 / \2 

KL> = ™l 2 +9l-M\ «> = ml + 9l-M 2 (15) 



and 



\m 2 Q x ) =rn 2 Q . + -g 2 B _ L (yzY 



(ml.) = m 2 Ui - -#|_ L (z/ 3 ) 2 , (m 2 d .) = m\ - -g 2 B ^ L (^) 2 (16) 

for i = 1,2,3. Note from the minus sign in the expressions for (m|.) and 
( m ui)y ( m di) that for this to be a stable vacuum the soft mass parameters 
m 2 L . and m 2 d . must always be positive at the B-L scale. This was shown 
to be the case. However, the same is not required for the m 2 2 , m 2 . and m^. 
parameters. These can become negative at the B-L scale as long as (m 2 ), 
(m^) and (mn ) are positive. This has important implications for bosonic 
superconductivity, as we will discuss in Section 5. For simplicity, we will 



always take m 2 to be positive at any scale, as was done in [30, 3 lj . 

Second, scale all parameters down to ~ 10 2 GeV. Here, one of the di- 
agonalized Higgs soft masses, indicated by a prime, becomes negative; that 
is, m 2 H , < 0. It follows that the up and down neutral Higgs fields develop 
non-vanishing VEVs given by 

^r-if?* ^ = ik {H ° } - (17) 

Evaluated at (H°), (H°) and (z/ 3 ), all other VEVs vanish. Therefore, this is 
a stable vacuum which, while continuing to break B-L symmetry at ~ TeV, 
now spontaneously breaks SU{2)l x U{1)y to U{1)em at the electroweak 
scale. This gives the Z and W vector bosons mass. Note that in our 
range of tan(3, (H°) (-^°)- Hence, although included in the numerical 
analysis, to simplify equations we will not display any (H°) contributions. 
For example, to leading order 

M z = V2{g 2 Y + giyi 2 (H«) . (18) 

The expressions for the slepton/squark masses at this minimum are each 
modified by additional terms proportional to the Higgs VEVs. By far the 
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largest such contribution is to the third family left- and right-handed up 
squark masses from their Yukawa interaction in (|6]). Ignoring much smaller 
D-term corrections, these are given by 



<«)) 



m 



1 2 



> 3 > 2 , («)> 



1 



111, 



where 



m 



U 3 ,u 3 



m: 



I I \ 



0\2 



•<(H») 



M 2 (19) 



(20) 



Comparing to (16), one sees that the mass parameters of U3 and 113 are 
modified by a positive Higgs VEV contribution. For this to be a stable 
vacuum, m 2 ^ must be positive at the electroweak scale. On the other hand, 
as long as ((mfj 3 )) is positive one can have m 2 ^ < 0. This is consistent with 
the conclusions at the B-L scale. However, to leading order 



: 1 2 
? 3 + t^9b-l 



(21) 



It follows that if ?71q 3 < 0, the potential is most destabilized in the D 3 direc- 
tion. Note that all other Higgs VEV contributions, either through F-terms or 
D-terms, are much smaller. Hence, with the exception of the splitting of the 
Us and D% mass parameters, all conclusions regarding soft masses reached 
at the B-L scale remain valid. For simplicity, we will no longer notation- 
ally distinguish between soft mass parameters m 2 and their Higgs corrected 
values m 2 . 

Finally, using the above results one can calculate the £>-L/electroweak 



hierarchy. It follows from (129) and (17) that 



(H°) ~ 9B-L 



Ki I 



(22) 



In the analysis of [221 EI], tan(3 was limited to 6.32 < tanf3 < 40 and there 
was a specific range of initial parameters. For a generic choice in this range, 
it was found that 



19.9 < 



fa) 
(F°> 



< 126 



(23) 



This demonstrates that a stable vacuum exists with an phenomenologically 
viable 5-L/electroweak hierarchy. Within this range of parameters, m 



< 



0. All other slepton/squark soft masses are positive with the exception of 
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7Bq 3 , which is negative. Hence, although the 5-L/electroweak vacuum is 
a minumum, at the origin of field space the potential is unstable in the D% 
direction. This has interesting applications to bosonic superconductivity and 
will be discussed in detail in Section 5. 

Recently, this analysis has been expanded to a much larger initial pa- 
rameter space using a completely numerical calculation of the RGE's. This 
will appear elsewhere. Suffice it here to say that, over this entire extended 
range, masses m^ 3 and m 2 H , are negative at the electroweak scale and induce 
a viable hierarchy of 

° m ~W)- ° m ' (24) 

However, within this expanded context, the squark/slepton masses are con- 
siderably less constrained. Specifically, each of the following combinations 
of soft scalar mass parameters at the electroweak scale can now occur: 1) 
all positive, 2) all positive except for m 2 ez < 0, 3) all positive except for 
m Q 3 < an d 4) combinations of these. We emphasize that in all cases the 
-B-L/electroweak vacuum is a stable absolute minimum of the potential and 
does not break color or charge symmetry. 



4 The B-L Cosmic String 

We begin by analyzing the theory at the B-L breaking scale. The preceding 
results show that this symmetry is radiatively broken by a VEV of the third 
right-handed sneutrino. Furthermore, evaluated at this vacuum, all squark, 
slepton and Higgs mass squares are positive. That is, this is a minimum of the 
potential energy and neither electroweak symmetry nor color is spontaneously 
broken at this scale. The situation at the origin of field space is more complex. 
As discussed above, it is possible for one or both of ra^ 3 and mg 3 to be 
negative. However, to introduce the basic cosmic string solution, in this 
section we analyze the theory assuming all soft mass parameters are positive. 
Under this assumption, the relevant physics is described by 

£ = \V^u 3 \ 2 - \f B - Lih ,F%_ l - V(u 3 ) , (25) 

where 

T> u .^ = - ig B ^ L A B - L ^ (26) 
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and 

v> 3 ) = <N 2 + %^N 4 . (27) 

The potential arises from two sources. The first term is the soft supersymme- 
try breaking third sneutrino mass term in (11 ) at the B-L scale. The second 
term arises as the pure third sneutrino part of the Db-l supersymmetric con- 
tribution in (|9]). Recall from the preceding RGE analysis that m 2 U3 = — |ra^J 
at the B-L scale. Hence, this potential is unstable at the origin and has a 
minimum at 

W 2 = • (28) 

9b-l 



Using this, potential (27) can be rewritten as 

^3) = %^(N 2 -<^3> 2 ) 2 . (29) 

Note that the soft supersymmetry breaking u 3 mass term has been re-expressed 
as the Fayet-Iliopoulos component of an effective D-term. It follows from this 
that the Higgs effect associated with (28) gives the Ab-l vector boson and 
the radial real scalar dv^ an identical mass 

m a b _ l = m Sv3 = V2g B -L(v?) ■ (30) 

The cosmic string solution to this theory is well-known [34J. Assuming 
a static solution that is translationally invariant in the z-coordinate, the 
cylindrically symmetric solution is of the form 

77 

^3 = e in > 3 )/(r) , A B - Lr = 0, A B -L* = a(r) . (31) 

9B-ir 

Here, integer n is the "winding number" of the string around the origin, 
which will always be assumed non-zero. The functions f(r) and a(r) have 
the boundary conditions 

/ ™ 1;/ ^ a ™ 1)C ^ (32) 

respectively. Before analyzing these functions further, it is important to note 
that there are two characteristic lengths associated with any cosmic string 
solution. These are 

r s = m£ , TV, = M- a \_ l . (33) 
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The explicit solutions for functions f(r) and a(r) will depend on the ratio 



r 2 



n = if ■ (3 4 ) 
In our case, we see from (30) that 



r. = r„ = (35) 

and, hence, 1Z is at the critical point 

K = 1 . (36) 

This is a consequence of the softly broken supersymmetry of our theory, and 
will have important implications when we study bosonic superconductivity. 
At the critical point, the equations for f(r) and a(r) simplify to 

/' = ^(l-a), - = n^) 2 (/ 2 -l) (37) 
r r \n\ 

where ' is the derivative with respect to r. Explicit solutions, even to these 
simplified equations, are not known, although their asymptotic expressions 
at small and large r have been evaluated [31]. However, precise numerical 
solutions for f(r) and a(r) exist in the literature, See, for example, [54]. 
We use both the asymptotic and numerical results throughout this paper. 
Another consequence of being at the critical point is that the energy density 
of the cosmic string simplifies to the exact result 

P = 2tt(z/ 3 ) 2 • (38) 

Let us now consider the theory at the electroweak breaking scale. As dis- 
cussed in the previous section, the up and down neutral Higgs fields develop 
non-vanishing VEVs given by 

<f/y = - (H°) = J-{H°) . (39) 
9y + 92 tan P 

Evaluated at (H°), (H°) and (z/ 3 ), all other VEVs vanish. Therefore, this is 
a stable vacuum which breaks both B-L and electroweak symmetries with a 
viable hierarchy. Does the electroweak phase transition effect the basic B-L 
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cosmic string solution? Since both Higgs field have vanishing B-L charge, z/ 3 
is electroweak neutral and the Yukawa coupling A„ 3 in ^ is of order 10~ 9 , 
the Higgs VEV contribution to the dynamical equations for z/3 is highly 
suppressed. Furthermore, it remains possible to choose initial parameters so 
that all soft masses are positive, even at the electroweak scale. It follows that 
the form of the cosmic string solution given in (31 ) does not change. Although 
the Higgs VEVs are no longer zero, since these fields are B-L neutral there are 
no topologically non-trivial solutions to the Higgs equations of motion [5T] . 
Henceforth, we assume that the Higgs fields are everywhere constants with 
the values given in ( 39 ) . 



5 Bosonic Superconductivity 

Bosonic superconductivity can occur if a charged scalar field develops a non- 
vanishing condensate in the core of the cosmic string [UJ . In the phenomeno- 
logical B-L MSSM theory discussed in this paper, there are a number of 
different ways that this could occur, each intricately related to other particle 
physics phenomena. Clearly, the first requirement for the existence of any 
such condensate is that a charged scalar mass squared at the origin of field 
space, that is, a soft supersymmetry breaking mass parameter plus small 
Higgs VEV corrections, becomes negative. As discussed in Section 3, there 
are several distinct ways in which this can occur. In this section, we exam- 
ine bosonic superconductivity within the core of the cosmic string in each of 
these scenarios. The entire analysis will be carried out at the electroweak 
scale. 

Case 1: All Soft Masses Positive 

This is the case described in the previous section. Since all soft supersym- 
metry breaking masses are positive, there can be no scalar condensates and, 
hence, no bosonic superconductivity at the core of the cosmic string. How- 
ever, such strings could exhibit fermionic superconductivity. This will be 
discussed in Section 6. 
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Case 2: Negative Soft Slepton Mass 

As discussed in Section 3, there is a region of initial parameter space such 
that, at the electroweak scale, all soft masses are positive with the exception 
of m 2 e3 < 0. This is the simplest case potentially admitting a non-zero conden- 
sate and, hence, we analyze it first. The relevant Lagrangian for discussing 
the vacuum of v 3 and e 3 is given by 

C = \T> V3 ^ 3 \ 2 - X -F B _ LiW F^ L + \V e3tl e 3 \ 2 - \F YfXU F^ - V(is 3 , e 3 ) (40) 
where 



V 



dp - igB-LAs-Ln, V e3ll = dfj, - ig B - L A B _ Lil - igvA Y ^ (41) 



and 



V(v 3 ,e 3 ) = m 



1 2 i 2 I |2 , 



9b-l 



(w 2 + \e3\r + 



9y 



e 3 



(42) 



The first two terms in the potential are the soft supersymmetry breaking 
mass terms in (11), while the third and fourth terms are supersymmetric 
and arise from the D B -l and Dy in ^ and ^ respectively. Contributions 
to (42) from the relevant Yukawa couplings in ^ are suppressed, since A„ 3 
and A es are of order 10~ 9 and 10 -2 respectively. Hence, we ignore them. The 
RG analysis tells us that both m 2 3 < 0, m 2 3 < at the electroweak scale. 
Hence, the potential is unstable at the origin of field space and has two other 
local extrema at 



rn 



1 ' 

9b-l 



(e 3 > = 



and 



(e 3 ) ; 



rn 



E3 



9b-l + 9y 



respectively. Using these, potential (|42|) can be rewritten as 

V(v 3 , e 3 



(43) 
(44) 



9b-l 



(N 



I \2\2 , 2 I 1 2 1 I 

w ) +g B -L\m [e 3 ] 



9b-l + 9y 



(\e 3 \ 2 - (e 3 )r ■ 



(45) 
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Let us analyze these two extrema. Both have positive masses in their 
radial directions. At the sneutrino vacuum (43), the mass squared in the 
direction is given by 



!<«*> 



9 2 B-M 2 -(9 2 B^L + 9Y)(e 3 y 



\m 



V3 I 



t-.i i 



(46) 



whereas at the stau vacuum (44), the mass squared in the z/3 direction is 



m 



<es> 



\m. 



S3 



2 (1 + ^)- 1 
9b-l 



\m 



v>, 1 



(47) 



Note that either (46) or (47) can be positive, but not both. To be consistent 
with the hierarchy solution, we want (43) to be a stable minimum. Hence, 

41 



we demand m 2 e3 | („ 3 ) > or, equivalently, that 



\m 



V3 1 



> \m 



es I 



(48) 



It follows from the RG analysis in [32] that one can always find a subregion of 
the initial parameter space so that this condition holds. We assume (48) for 
the remainder of this subsection. It then follows from (47) that m^ 3 |( e3 ) < 



and, hence, the stau extremum (|44|) is a saddle point. As a consistency check, 

(49) 



note that V\^ U3 ) < if and only if 

or, equivalently, 



V3 1 



> |m e3 r(l + -p— ) 



9y n-1/2 



9 2 b-l 



(50) 



This follows immediately from constraint (48). Finally, note that the poten- 
tial descends monotonically along a path C from the saddle point at (44) to 

dV 



the absolute minimum at (43). Solving the d 
found to be 



equation, this curve is 



Me = ((e 3 } ; 



N 2 (i + 



9y \-l\l/2 



9 2 B-L 



(51) 



Note that it begins at (63) for z/ 3 = and continues until it tangentially 



intersects the = axis at I1/30I 



l m e 3 l 

\m v A 



(z/3). From here, the path continues 



down this axis to the stable minimum at (43). 

We conclude that at the electroweak scale the absolute minimum of poten- 
tial (42) occurs at the sneutrino vacuum given in (43). The sneutrino scalar 
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can develop a non-zero winding around some point in three-space, leading 
to a cosmic string. Away from the core, this will still be described by the 
simple cosmic string solution in the previous section. Recall, however, that 
non-zero winding forces the function f(r) and, hence, u 3 to vanish at r = 0. 
This was not an issue for the simple cosmic string, since it was assumed that 
all squark/slepton masses were positive at the origin of field space. In the 
present scenario, however, the mass squared of e 3 , 

m l 3 \v 3 = m % + g 2 B-L"3 2 , (52) 

becomes negative as u s approaches the origin of field space. This potentially 
destabilizes the e 3 field in the core of the string, producing a scalar conden- 
sate. Whether or not this can occur is dependent on the relative magnitudes 
of the spatial gradient and the potential energy, which tend to stabilize and 
destabilize e 3 respectively. To analyze this, one can look at small fluctuations 
of e 3 around zero in the background of the simple z/ 3 cosmic string solution 
in Section 4. The equation of motion for e 3 is given, to linear order, by 

(d»d» + 2ig B 

Hol-M 2 - (9%-l + ^)(e 3 ) 2 )e 3 = , (53) 



where v% and Ab-l^ were defined in (31). Using the ansatz 

e 3 = e-*e 30 (r) , (54) 



equation (53) simplifies to 

&>_ _ ld_ 
dr 2 r dr 



e 30 + ^ e so = w e 30 (55) 



where 



nr)=gl^) 2 f(rf-( 9 l„ L + gl)(e 3 )> + n^ . (56) 



Note that we have chosen e 30 in (54) to be a function of radial coordinate r 
only and, hence, not to wind around the origin. If this two-dimensional 
Sturm-Liouville equation has at least one negative eigenvalue, the corre- 
sponding u becomes imaginary. This destabilizes e 3 , implying the existence 
of an e3 condensate in the core of the cosmic string. 
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Case 3: Negative Soft Squark Mass 

As discussed in Section 3, there is a region of initial parameter space such 
that, at the electroweak scale, all soft masses are positive with the exception 
of rriQ 3 < 0. The electroweak phase transition breaks the left-handed SU(2) L 
doublet Q 3 into its up- and down- quark components U 3 and D 3 respectively. 
The leading order contribution of the Higgs VEVs to their mass splits the 
degeneracy between these two fields, destabilizing the potential most strongly 
in the D 3 direction. For this reason, the relevant Lagrangian for analyzing 
this vacuum can be restricted to 



£ - \V U3tl is 3 \ 2 - -Fb-l^vFb-l + \V D ^D 3 \ 2 - -Fy^Fy 



where 



~^ F su(2)^uF^ {2) - ^F S u(3)^uF^ {3) - V(v 3 , D 3 ) (57) 



= 0^- igB-LA B ~Lti , (58) 

^. r, -9B-L a -9Y A -A -A 



and 



V(v 3 ,D 3 ) = ml,\v 3 \ 2 + m 2 Dz \D 3 \ 2 + ^(N 2 + l\D 3 \ 2 ) 2 (59) 

+ i + i)|£>, 

2 V 36 4 3 n 1 



,2 

~2 ' " V 

i i 



The first two terms in the potential are the soft supersymmetry breaking 



mass terms in (11), while the remaining terms are supersymmetric and arise 
from the D B -l, D y in @,([8} and D S u(2) L , D S u(3) c respectively. Using 
Ad 3 ~ 5 x 10~ 2 , the hierarchy given in (24) and assuming \mjj 3 \ is of order 
|m^ 3 |, terms proportional to the Higgs VEVs are small and are ignored in 
(59). For simplicity, we henceforth drop the small g B _ L /9+g Y /36 piece of the 
D-term contribution. The RG analysis tells us that both m 2 3 < 0, mf^ < 
at the electroweak scale. Hence, the potential is unstable at the origin of 
field space and has two other local extrema at 



9 

m. 



W 2 = --^S (D 3 ) = 0, (60) 
9b-l 
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and 



(1/3) = 0, (D 3 y 



m 



D 3 



respectively. Using these, potential (]59|) can be rewritten as 

V(u 3 ,D 3 



9b-l 



+ ft a /4 + g l/3 , 



9b-l , 
3 1 

(A.) 3 ) 2 



(61) 



(62) 



Let us analyze these two extrema. Both have positive masses in their 
radial directions. At the sneutrino vacuum (60), the mass squared in the D 3 
direction is given by 



m D :i \(u 3 ) 



9 2 B-L 



■9l , gl 



\m D3 



(63) 



whereas at the stau vacuum (61), the mass squared in the z/ 3 direction is 



9 2 B-L 



(D3) 2 - g 2 B- L {^) 2 = \ m D s 



9 2 B-L 



\m 



^3 1 



(64) 



Note that either (63) or (64) can be positive, but not both. To be consistent 
with the hierarchy solution, we want (60) to be a stable minimum. Hence, 
we demand tn 2 D3 \^ 3 ) > or, equivalently, that 



\m 



V3 I 



> 3\mD 3 



(65) 



The RG analysis in [31] shows that one can always find a region of the 
initial parameter space so that this condition holds. We assume (65) for the 
remainder of this subsection. It then follows from (64) that m 



(D 3 ) 



< and, 

hence, the stau extremum (61 ) is a saddle point. As a consistency check, note 
that V\/ U3 \ < ^|(d 3 ) if and only if 



g 2 B -M A >( g j + 9 j)(DsY 



or, equivalently, 



\m Vz 1 > \m D3 



g 2 B-L 



1/2 



(66) 



(67) 
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This follows immediately from constraint (65). 



We conclude that at the electroweak scale the absolute minimum of po- 
tential (59) occurs at the sneutrino vacuum given in (60). The sneutrino 



scalar can develop a non-zero winding around some point in three-space, 
leading to a cosmic string. This is described, away from the core, by the 
simple cosmic string solution in the previous section. Recall, however, that 
non-zero winding forces the the function f(r) and, hence, v% to vanish at 
r = 0. In the present scenario, the mass squared of D%, 



m 



m 



D 3 



+ 



9%- 



L 2 

-fa 



(68) 



becomes negative as u 3 approaches the origin of field space. This potentially 
destabilizes the D3 field in the core of the string, producing a scalar conden- 
sate. To analyze this, one can look at small fluctuations of D3 around zero 
in the background of the simple v% cosmic string solution. The equation of 
motion for D 3 is given, to linear order, by 



(d^ + 2i^A B _ LO d e 



9b-l 
9 



A 



A M 

B-L/i A B-L 



1 (9b-l 1 12 
+(— s— "3 



;| + |)(D 3 m = o 



where v 3 and Ab-l^ were defined in (31). Using the ansatz 



D 3 = e^Dso(r) 



equation (69) simplifies to 

_ ld_ 
dr 2 r dr 



)D 30 + VD30 = w 2 D 30 



where 



V-(r) = %^3> 2 /(r) 2 - (f + |)(^) 2 + 



n 2 a{r) 2 
~9 r 2 



(69) 



(70) 



(71) 



(72) 



Note that we have chosen D30 in ( 70 ) to be a function of radial coordinate r 
only and, hence, not to wind around the origin. If this two-dimensional 
Sturm-Liouville equation has at least one negative eigenvalue, the corre- 
sponding u becomes imaginary. This destabilizes D 3 , implying the existence 
of an D3 condensate in the core of the cosmic string. 
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Numerical Analysis of Boson Condensates 

Let us analyze the formation of a scalar condensate in a more general setting. 
Consider a U(l) x U(l) gauge theory with two complex scalar fields and 
a charged under the gauge group as q^ = 0, q^ 7^ and q a ^ 0, q a ^ 
respectively. U{1) and U(l) are motivated by Uy and Ub-l m the previous 
sections. Similarly, scalar corresponds to the right-handed sneutrino z/3. 
A condensate can potentially form in the a field. Unlike previous analyses 
in the literature, here, in addition to the usual U(l) charge of a, q a is also 
non- vanishing. This is motivated by the fact that all squarks and sleptons 
in the B-L MSSM theory carry non- vanishing B-L charge. After finding the 
necessary conditions for a condensate to form, we will apply the results to 
the specific cases discussed above. 

The Lagrangian density for this generic theory is given by 

C = \V^\ 2 - \f,^ u + |X>| 2 - - V(<P, a) (73) 



where 
and 



i>n = dp - iq^gA^, V^ = d^- iq^gA^ - iq a gA^ (74) 

V(<f>, a) = ^(|0| 2 - t/J) 2 + m V| 2 + ^(|a| 2 - t/ 2 ) 2 . (75) 

The coefficients \^,X a and (3 are chosen to be positive. Potential (75) has an 
extremum at (0) = (a) = 0. If one chooses the coefficients so that the 
effective a mass squared at this extremum is positive, that is, 



then (0) = 7/0, (a) = is a local minimum. This vacuum spontaneously 
breaks the U(l) symmetry and admits a cosmic string solution in of 



the form discussed in Section 4. Potential (75) has a second extremum at 
(0) = 0, (a) = f]cr. This may or may not be a local minimum of the poten- 
tial depending on the choice of parameters. In all cases, however, one can 
constrain the cosmic string vacuum to be deeper than the a extremum by 
choosing 

V7$ > hvt , (77) 

which we do henceforth. 
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As discussed in Section 4, there is a U(l) cosmic string solution of the 
associated <fi and equations of motion given by 



And 



V<f>f( r ) 



0, A 6 



n 

air) 



(78) 



where integer n is the "winding number" of the string around the origin. The 



functions f(r) and a(r) have the boundary conditions given in (32). In the 
theory we are considering, the effective mass squared of a at arbitrary is 



2 



(79) 



This becomes negative as approaches the origin of field space, potentially 
destabilizing the a field in the core of the string and producing a scalar 
condensate. Such a condensate would break both U(l) and U(l) symmetry. 
Whether or not this can occur is dependent on the relative magnitudes of 
the spatial gradient and the potential energy, which tend to stabilize and 
destabilize a respectively. To analyze this, one can look at small fluctuations 
of a around zero in the background of the simple cosmic string solution. 
The equation of motion for a is given, to linear order, by 



(dpdP + 2iq a gA e d e - &A^)a 



+m\ 2 - 



a 



(80) 



where and A a were defined in ( 78 ) . Using the ansatz 



equation (80) simplifies to 



d 2 



1 d 



where 



V(r) = Pvlf(r) 



I (T + V(T 



wV 



2 at r 2 



(81) 



(82) 



(83) 



Note that we have chosen ctq in (81) to be a function of radial coordinate r 



only and, hence, not to wind around the origin. We want to emphasize the 
term in (83) proportional to a 2 /r 2 . This appears precisely because the a field 
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has non-vanishing charge under U(l) as well as under U(l), a situation not 
previously discussed in the literature. However, for the reasons mentioned 
above, it must be included in the analysis of this paper. 

If this two-dimensional Sturm-Liouville equation has at least one negative 
eigenvalue, the corresponding u becomes imaginary. This destabilizes a, 
implying the existence of a a condensate in the core of the cosmic string. 
Note that if the condensate was not charged under U(l), then the last term 



in (83) would not appear and potential V(r) would be monotonic. As was 
discussed in jl?], for sufficiently small m 2 ^^ a condensate will always form 
under these conditions. However, a is charged under U{1) and, hence, the 
a 2 /r 2 term in (83) must be included. Since this term is positive and provides 



a repellent force for large r, it may prevent a bound state from forming. 
For the remainder of this section, we will discuss the results of a numerical 



solution to the Sturm-Liouville equation (82) with potential (83). The details 
of this solution are presented in the Appendix. 

The a 2 /r 2 term is smallest and, hence, the least disruptive to the forma- 
tion of a condensate for winding number n — 1. Therefore, we carry out the 
analysis for the singly wound cosmic string. Motivated by the softly broken 
supersymmetric B-L MSSM theory, the calculation will be further restricted 
in two ways. First, take q 2 = q\. Second, we constrain the parameters to the 
critical coupling point where 



A 
2g 



* l - ( 84 ) 



It follows that the functions f(r) and a{r) simplify to solutions of (37). These 
equations can be solved numerically [M] , and we input them into our analysis 
of the Sturm-Liouville equation. Finally, simplification can be achieved if we 
take 

= (85) 



thus setting m 2 at the cosmic string vacuum, given in (76), to zero. Potential 
V(r) then becomes 

V(r) = /3^(/(r) 2 - 1) + ■ (86) 

Should a negative energy bound state exist for some choice of /3, the conden- 
sate will persist if we continuously deform m 2 ]^^ away from zero to a small 
positive value. 
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In units where is one, the Sturm-Liouville equation (82) with potential 



(86 ) depends on the single parameter (3. Explicit solutions of this equation for 
several values of f3 are shown in Figure 1. For each of these values, a negative 
energy eigenvalue and normalizable bound state wavefunction exists and are 
shown in the Figure. As discussed in the Appendix, we find that a negative 
energy eigenvalue will exist for any 

(3 > (3 critical ~ 0.42 . (87) 



Hence, for sufficiently large (3 satisfying (87), that is, for sufficiently deep 
potential, the a field is destabilized and a non-vanishing a condensate will 
form. However, for (3 less than Perineal the eigenvalue becomes positive and 
the wavefunction oscillatory, signaling a meta-stable solution. Hence, for 
(3 < 0.42 the potential is not sufficiently deep and a a condensate will not 
form. These results can immediately be applied to the B-L MSSM theory 
with a negative soft right-handed slepton mass m 2 < described in Case 2 



above. Identifying <p = z/ 3 , a = e 3 and comparing (74), (75) to (41), (45) using 



(g3j,@, we find that 

P = 9b-l, ^ = 2 9b-l, ^ = 2 (9b-l + 9y)> 

(88) 

In particular, evaluated at the electroweak scale 

(3 = g\_ L ~ 0.1075 < 0.42 (89) 

suggesting the absense of a bosonic condensate in the core of the cosmic 
string. 

The formalism applicable to the B-L MSSM theory with a negative soft 
left-handed squark mass squared requires a change in the relative charges 
of and a. Instead of taking q 2 = cfy as we did previously, now choose 
q\ = 9g^. For n = 1 winding number at the critical point and vanishing 



m^L., the Sturm-Liouville equation determining the a condensate is (82) 



with potential (83) now given by 



V{r) = /3^(/(r) 2 - 1) + . (90) 

A numerical analysis completely analogous to the one used above, leads to 
the conclusion that a non- vanishing a condensate will occur for 

/? > Pcritical ^ 0.14 . (91) 
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f3 = 0.5 



ul = -0.0231 



Figure 1: Negative energy ground state solutions of the slepton stability 
equation for (3 = 0.8 and (3 = 0.5 respectively. The energy eigenvalues cu 2 . 
are shown as red lines, with the associated normalizable wave functions o"o 
depicted in blue. Note the positive "bump" in potential V due to the a 2 /r 2 
term. 



Note that this is smaller than the previous bound given in (87). This is 



consistent with expectations since the destabilizing a 2 /r 2 term is now smaller 
by a factor of 9. For (3 less than ftcriticah however, the eigenvalue becomes 
positive and the wavefunction oscillatory, signaling a meta-stable solution. 
Hence, for (3 < 0.14 the potential is not sufficiently deep and a a condensate 
will not form. By construction, these results can immediately be applied to 
the B-L MSSM theory with a negative soft left-handed squark mass m|, 3 < 
described in Case 3 above. Identifying = z/ 3 , a = D 3 and comparing 



(74), (75) to (58), (62) using (60), (61), we find that 







9 2 b-l 



A, 



m. 



2<7 



2 ' 

9b-l 



B-Li 



4 3 



ni 



D 3 



S 2 74 + gi/3 



(92) 
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In particular, evaluated at the electroweak scale 

(3 = 9 -Z-k ~ 0.0358 < 0.14 (93) 

suggesting the absense of a bosonic D 3 condensate in the core of the cosmic 
string. 



6 Fermionic Superconductivity 

A second possible signature of cosmic strings is superconductivity arising, not 
from boson condensates but, rather, from zero-modes of charged fermions |47J. 
The relevant fermions are those with bilinear couplings to a scalar field that 1) 
has a non-vanishing VEV at radial infinity and 2) which winds non-trivially 
around the center of the string. Exactly which fermions, if any, develop zero- 
modes is dependent on the theory under consideration and on the explicit 
cosmic string background (531 S3 ED]- In the B-L MSSM theory described 
in this paper, the structure and properties of potential zero-modes is very 
specific. 

The cosmic string background described in Section 4 is constructed from 
the right-handed sneutrino, which has non- vanished VEV at radial infinity 
and non- zero winding n around the origin. We begin, therefore, by consider- 
ing the fermions which couple to it. In the B-L MSSM theory, the coupling 
of the right-handed sneutrino u 3 to charged fermions is completely specified 
by the superpotential 

W = --- + K 3 L 3 Hv 3 , (94) 

where X U3 is the third family neutrino Yukawa coupling of order 10~ 9 . It 
follows that the relevant physics is described by 



-K z {^ LE -i ) LH+^ + hc) + ... , (95) 

where = + igB-L^B-Ly,- The associated equations of motion are 

ia^V Ell ip LE - - K^lh+vI = 0, 
ia^d^ LH + - X U3 ip LE -u 3 = 0. (96) 

We want to solve these in the background of the cosmic string defined by the 



transverse functions z/3 and Ab-l^ in (31 ). Therefore, first consider solutions 
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of (96) that are independent of time and the ^-coordinate. Denoting these 



transverse fermions by /3 LE -(x,y) and ^LH+(x,y), equations (96) become 

id- l T>i/3 LE - - X^Plh+Vs = °> 
iv%P L H+ - X„ 3 P LE -v 3 = (97) 

where V Efl = 9^ + ^b-lAb-l/i and % — 1,2. It follows from the index of 



the corresponding Dirac operator that (97) has \n\ linearly independent pairs 
of normalizable zero-modes jlSl ED]. These modes are eigenstates of the 03 
operator, 

a 3 (3 = . (98) 
\n\ 

Thus, for a given n the zero-modes have the same chirality and each is de- 
scribed by a complex scalar function. Specifically, in cylindrical coordinates 
the solutions are [IHl [50] 

faH+(r,0) = U l LHt (r)e« l+ ^ 9 (99) 

where — f+|<£<f— §. The radial functions can be explicitly evaluated 
asymptotically. As r — > 0, one finds 

U l LE -{r) ~ r-'+t-i, 

~ r i+ t~^ (100) 

whereas for r — > 00 

^W (r) ' ^ 3+ W~e-^<-3>r. (1Q1) 
Due to the exponential decay, the range of the fermionic solutions is of order 

r <~Kk)- (102) 



Note from (35) that for our specific theory the radius of the cosmic string is 



r„ ~ 7— r and, hence, 



r f - 9B-L „ 1q9 (1Q3) 



r X 

1 S /y U3 
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That is, the radial zero-mode solutions are 10 9 times wider than the vortex 
core. This indicates the extremely diffuse nature of the fermionic solutions 
in the B-L MSSM theory. Any normalizable solution with general boundary 
conditions is a linear combination of these zero-modes. We refer the reader 
to [1H1 EDI E3] for a detailed derivation of these properties. 



Now consider full four- dimensional solutions of (96) of the form 



^LE- = 0LE-( X ^y) a ( Z ^) , ( 104 ) 

4>lh+ = /3 LH +{x,y)a(z,t)* . (105) 



Since (3 LE -(x,y) and /3lh+(x,v) satisfy the transverse Dirac equations (97) 



it follows from (96) and (98) that 

/ n d \ 

G»-H&H 2 ' i)=0 - (106) 

Without loss of generality, we henceforth assume that the winding number is 



positive. Then (106) implies that a(z,t) = f(z + t). Thus, an n > winding 



of the z/3 scalar in the B-L cosmic string solution (31) induces left-moving 
fermionic currents in the cosmic string composed of i>LE~ an d 4>lh+- 

Anomaly cancellation on the string worldsheet [UJ requires that there 
be chiral fermions which couple to a scalar which winds oppositely to z/3. 
Since the theory is supersymmetric, this cannot be the conjugate field z/3. 
Furthermore, since the Higgs fields are neutral under B-L transformations, 
they cannot have topologically stable winding around the string core even 
though they have non- vanishing VEVs. Note, however, that if the left-handed 
sneutrino Ne 3 gets an expectation value, then it follows from the equations 
of motion that this must wind oppositely to z/3. Any chiral fermions coupling 
to the wound solution Ne 3 will then generate right-moving currents on the 
cosmic string, canceling the anomaly. Does A3 have a non-zero expectation 
value? The answer is affirmative, as we now show. 

The neutral scalar fields in the B-L MSSM theory are the H°, H° com- 
ponents of the Higgs fields and the left- and right-handed sneutrinos Aj, V{ 
for % = 1, 2, 3. Since z/3 is the only right-handed sneutrino to get a non- zero 
expectation value, we need only consider the third family. The potential 
energy of these neutral fields is found to be 

V = V F + Vb—L +Vy + V SU (2) + Vsoft (107) 
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where 

V F = J2\Fm\ 2 = A^(|z, 3 ri^T + l^| 2 |iV3r + |iV3| 2 |if | 2 ) (108) 

m 

+ /i 2 (|#°| 2 + \H°\ 2 ) - \u 3 (^3N 3 H° + he) , 
Vb-l = \d\_ l = %^((N 2 - |iV 3 | 2 ) 2 , (109) 

v Y = \d\ = 9 f(\Hr - \sy - , (no) 

2 



l#°l 2 - 


\N 3 \ 


\H°\ 2 + 


\H°\ 



Vsum = \ D2 su { 2) = f H#T + m 2 + l^l 2 ) 2 , (m) 

Vsoft = mljN 3 \ 2 + ml 3 \u 3 \ 2 + m 2 H \H \ 2 + ml\H \ 2 (112) 
- (BH°H° + he) + (A U3 u 3 N 3 H° + he) . 

Let us solve for the expectation values for each neutral scalar subject to the 
hierarchy condition 

(Ar 3 ><<(tf >,(tf°)<< W ■ (H3) 
The dV /dv 3 = and dV /dH° = 0, dV /dH° = equations lead to the 



non-zero expectation values for (v 3 ) and (H ) , (H ) presented in (129) and 



(17) respectively. The dVo/dN 3 = equation then gives 



,, M (KttfF) - <ui\ 

{N3) - mh-fi-M* ■ (U4) 

Therefore, following the electroweak phase transition the left-handed sneu- 
trino acquires a very small expectation value. For example, assuming ~ 
(H°), A V3 ~ \ V3 (H°) and m% a - g 2 B - L (v 3 ) 2 ~ g 2 B _ L {u 3 ) 2 , it follows that 

{iV 3 > ~ (10- 10 - 10- 12 )(^ 3 > (115) 



for the i?-L/electroweak hierarchy given in (24). That is, (N 3 ) is on the 
order of the neutrino masses. This is sufficient, however, to provide the 
right-moving fermionic zero-modes on the cosmic string required by anomaly 
cancellation. Note that the vanishing of (z/ 3 ) at the center of the cosmic 
string will set (N 3 ) = 0, consistent with a non-trivial winding of N3. 



27 



To see how these arise, note that the coupling of the left-handed sneutrino 
N% to charged fermions is specified by the superpotential 



W 



+ X e3 L 3 He 



3 • 



(116) 



where A e3 is the third family r Yukawa coupling of order 5 x 10 2 . The 
relevant physics is then described by 



£ = i ^Le+ aflT> e^ Le + + ilpLH-^d^LH- 
-^eMbe+^LH-Nz + he) + . . . 



117) 



where V efJL = <9 M — igB-L^B-L^- It follows from the above analysis that with 



respect to the cosmic string defined by u 3 , A b _l/i in (31) and the associ- 
ated background N 3 , there will be \n\ linearly independent pairs (3 Le +(x,y), 
0lh- {x, y) of transverse normalizable fermion zero-modes. These modes are 



eigenstates of a 3 and have a structure similar to (99). However, whereas the 



solutions in (99) and their o~ 3 eigenvalue are indexed by n, the winding of 
z/3, these zero-modes are indexed by —n, the winding of N3. The small r 



behaviour remains similar to that in (100). Now, however, as r — > 00 



~ e 



-Aeo<A r 3>r' 



and, hence, the range of these fermionic solutions is of order 

1 

r F 



It follows from this, ( |T02[ ) and ( pL15[ ) that 

r F ~ (10 3 - 10 5 )r/ . 



(118) 



(119) 



(120) 



Therefore, these radial zero-mode solutions are even more diffuse around the 
cosmic string core. The full four-dimensional solutions are again of the form 
i>L = (3ij(x,y)(y(z,t). Now, however, the function a satisfies 



Of 



d_ 

dzJ 



)a(z,t) 







(121) 



implying that a(z, t) = f(z — t). Thus, the n < winding of the N 3 scalar in 



the background of (31 ) induces right-moving fermionic currents in the cosmic 
string composed of ?p Le + and ipiH- ■ 
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Having found both the left- and right-moving charged fermionic modes of 
a B-L MSSM cosmic string, we want to analyze whether they can lead to cos- 
mologically observable phenomenon and, specifically, to superconductivity. 
As can be seen from the above discussion, the electroweak phase transition is 
necessary to produce the right-moving modes to cancel the anomaly How- 
ever, the superpotential term in (116) will then generate, in addition to the 
fermion coupling to N 3 discussed above, a Yukawa mass \ e3 ip LE -{H )%lj Le + 
for the tauon. In addition, a non-zero /x-term /iHH is required in the super- 
potential to make Higgsinos massive. Specifically, a mass term of the form 
HipH+ipH- will appear. Both masses are much larger than neutrino masses 
and, as a result, the zero-modes of the previous discussion will be lifted; 
generically, turning into massive bound states jlHl SHI E2 E2] with mass 
m T = 1.776 GeV and mmggsino ~ H respectively. For an applied electric field 
in the string frame satisfying E <C 27rm 2 , massive bound states cannot form 
persistent currents. Their maximum electric current is given by 



Note that J ma x is directly proportional to the applied field. As soon as E = 0, 
this charged fermionic current will relax to zero. The maximal such currents 
generated in a B-L MSSM cosmic string by tauon bound states, for example, 
would be of order 10 6 -10 7 A For E 3> 27rm 2 , it is possible to obtain currents 
in this setting which are close to superconducting. However, for that to 
happen one needs to be in the regime 

B(v/c) > 10 3 (m/leV) 2 Gauss , (123) 

where v is the transverse string velocity. Hence, for a superconducting current 
of tauons, the required magnetic field would be at least of order 10 21 Gauss, 
far larger than any observed cosmological B field. Thus, fermionic currents of 
the B-L MSSM cosmic string are unlikely to have any observable cosmological 
effects. 
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Appendix: Numerical Analysis of the Stability 
Equation 

In this Appendix, we present a numerical procedure for determining the exis- 
tence of negative eigenvalue, normalizable solutions of the stability equation 



d 2 



1 d 



dr 2 r dr 



i o- + Vvo 



u 2 a 



where 



q 2 a{r) 2 



(124) 



(125) 



Although applicable in general, we specify our algorithm for the simplest 
case discussed in the text where n = 1, q 2 = qi and m 2 a = 0, resulting in 



(85). Potential (125) then simplifies to 



V(r) = P(f(rf 



D + 



a(r)' 



(126) 



where we have set r]^ = 1 and, hence, the radial coordinate r and parameter 
(3 are dimensionless. Note that V depends only on the single parameter (3. 
Furthermore, impose the critical coupling constraint ^-p = 1, thus simplifying 
the functions f(r) and a(r) to be solutions of (37). These equations have 



been solved numerically in the literature [51] and we input their solutions 
directly into our analysis of the stability equation. 

To prove the existence of a boson condensate for a fixed value of parameter 



(3 in (126), it suffices to find a negative energy ground state solution to (124) 



Hence, one can impose the boundary conditions 



00 



r=0 



d r a 



. 



r=0 



In addition, to ensure that <J (r) is normalizable constrain 



. 



(127) 



'128) 



Note that if a bound state exists, its eigenvalue can never be more negative 
than the depth of the potential energy. Hence, the possible range of values 
for bj 2 is limited to 



-/3 < u 2 < 0. 



(129) 
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For each such u 2 , there is a solution ctq which satisfies (124) with boundary 



conditions (127). Generically, however, this solution will not be normalizable. 



To see this, note that for large r the solution to (124) must be of the form 



-|w|r 



C 2 r 



-l/2 e |w|r 



(130) 



where C\ and C 2 are continuous functions of u 2 . If for a chosen value of u 2 
C 2 is non-vanishing, then o"o diverges at large r and the renormabiliability 
constraint (128) is not satisfied. Note that C 2 ^ can be either positive of 
negative. If positive, o"o r —^ +oo, that is, the wavefunction "flips up" at large 
r. On the other hand, if C 2 is negative, (Tq — > — oo and the wavefunction 
"flips down". It is only if C 2 exactly vanishes for a specific value u 2 in (129), 



that constraint (128) is satisfied and the wavefunction normalizable. 

Two scenarios are then possible. First, if when u 2 is varied over the 



entire range (129) C 2 is always greater than, or always less than, zero, then 



the wavefunction is never normalizable and a negative eigenvalue ground 
state does not exist for this choice of parameter /3. Second, if when u 2 is 



varied over range ( 129 ) C 2 changes sign, then there must be an u 2 for which 



(131) 



CK 2 ) = o , 

since C 2 is a continuous function of u 2 . Hence, for this choice of (3 a normaliz- 
able ground state solution for <To exists with negative energy Wq. These results 
give us an explicit algorithm for computing the existence, or non-existence, 
of a boson condensate. This is: 



Choose a fixed value for parameter /?. 
Vary cu 2 over the range (129). 



For each value of lo , numerically solve (124), (126) for the ground state 



wavefunction cr satisfying boundary conditions (127). We do this by 



implementing the Runge-Kutta method on Mathematica. 
Plot do versus r for all values of u 2 . 

If all these curves "flip up" or "flip down", then there is no negative 
energy ground state. However, if these curves "flip up" for small val- 
ues of bj 2 but "flip down" for larger values, then a ground state with 
negative energy ujq does exist. 
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6. To compute and the associated normalizable wavef unction, we nu- 
merically identify the interval which contains u 2 . We then iterate this 
procedure until we obtain the ground state energy and wavefunction to 
the desired precision, thus approximating the solution to the stability 
equation. 



To make this concrete, in |Figure 2] we carry out this algorithm explicitly 
for parameter (3 = 0.8. Observe that o"o "flips up" for small u 2 , but "flips 
down" for larger values of u 2 . This signals the existence of a negative energy 
ground state occurring in between, when <jq — — £ 0. The numerical value of 
u 2 = -0.1421 and 
Figure. Note that 



u 2 = —0.1421 and the normalizable wavefunction are both indicated in the 



= 0.1776 , (132) 



that is, the bound state energy is 17.76 % of the depth of the potential. 
Let us now carry out this computation for smaller values of (3. The results 



for (3 = 0.5 are shown in Figure 3 Again, note that o"o "flips up" for small 
u 2 , but "flips down" for larger values of u 2 . This signals the existence of 
a negative energy ground state occurring in between, when <7o — — £ 0. The 
numerical value of u 2 = —0.0231 and the normalizable wavefunction are both 
indicated in the Figure. In this case, 

= 0.0462 , (133) 

that is, the bound state energy is 4.62 % of the depth of the potential. Note 
that the percentage size of the eigenvalue relative to the depth of the potential 
has substantially decreased over the = 0.8 case above. This indicates that 
for some value of f3 not too much smaller than 0.5 a negative energy ground 
state might cease to exist. To explore this further, we apply our algorithm 
to a range of values of parameter (3. The ground state energy for each (3, 
as well as their fractional depth with respect to the potential, are shown 
in Table 1 Note that as (3 approaches ~ 0.42, uj 2 — > and is a rapidly 



decreasing percentage of the potential depth. Indeed, we find that for 

(3 < (3 critical ~ 0.42 , (134) 



there is no negative energy bound state solution to (124),(126). Two concrete 



examples of this are (3 = 0.35 and /3 = 0.1. Our numerical results for these 
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l^ol 2 




1 


0.2404 


0.2404 


0.8 


0.1421 


0.1776 


0.7 


0.0977 


0.1396 


0.5 


0.0231 


0.0462 


0.45 


0.0094 


0.0209 


0.42 


0.0027 


0.0064 



Table 1: The ground state energy corresponding to different values of j3. 
Note that as the potential becomes more shallow, the ground state energy 
decreases relative to the depth of the potential. 



parameters are shown in Figure 4 and Figure 5 respectively. For both cases 
we see that, unlike the previous examples, a always "flips up " for all values 



of a; satisfying (129). It follows that in each case there is no negative energy 



ground state solution to the stability equation. To conclude: we have shown 



numerically that the stability equation (124) with potential (126) admits a 



negative energy ground state normalizable solution if and only if 



0># 



critical 



0.42 



(135) 



As discussed in the text, a similar analysis must be carried out with the 
charges chosen to be eft 



9q*. This changes potential (126) to 



V(r) = (3(f(r) 



D + 



9r 2 



(136) 



The numerical analysis of this case gives the same qualitative results, so we 
won't present it here. Suffice it to say that, due to the weaker repulsion term 
in the potential, the critical value for (3 is lowered. Specifically, we find that 



the stability equation (124) with potential (136) will admit a negative energy 



normalizable ground state if and only if 



P>P, 



critical 



0.14 



(137) 
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Figure 2: Family of <To solutions for the initial value problem with j3 = 0.8 
and u 2 varying from -0.12 to -0.156. Note that the asymptotic behaviour of 
the wavefunction changes sign. The ground state occurs at u 2 = —0.1421 
and its associated normalizable ground state is indicated in red. 




ti) 2 = -0.036 



Figure 3: Family of a® solutions for the initial value problem with = 0.5 and 
u 2 varying from -0.003 to -0.036. Note the changing sign in the asymptotic 
behaviour of the wavefunction. The ground state occurs at u 2 = —0.0231 
and the associated normalizable ground state is indicated in red. 
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/3=0.35 



5 10 15 20 

Figure 4: Family of <j solutions for the initial value problem with = 0.35 
over the entire allowed range of u 2 . Note that the asymptotic values of the 
wavefunctions are always positive, diverging to +oo. This corresponds to the 
stability equation admitting no negative energy ground state. 




0.75 - 
0.5 

£=0.1 

0.25 : 



2 4 6 8 10 

Figure 5: Family of do solutions for the initial value problem with = 0.1 
over the entire allowed range of u 2 . Note that the asymptotic values of the 
wavefunctions are always positive, diverging to +oo. This corresponds to the 
stability equation admitting no negative energy ground state. 
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